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Abstract. We study the spectral properties of the transition semigroup of the killed one- 
dimensional Cauchy process on the half-line (0,oo) and the interval (—1, 1). This process 



is related to the square root of one-dimensional Laplacian A = with a Dirichlet 

exterior condition (on a complement of a domain), and to a mixed Steklov problem in the 
half-plane. For the half-line, an explicit formula for generalized eigenfunctions ?p\ of A is 
derived, and then used to construct spectral representation of A. Explicit formulas for the 
transition density of the killed Cauchy process in the half-line (or the heat kernel of A in 
(0, oo)), and for the distribution of the first exit time from the half-line follow. The formula 
for ip\ is also used to construct approximations to eigenfunctions of A in the interval. For 
the eigenvalues A„ of A in the interval the asymptotic formula X n = ^ — ) is 

derived, and all eigenvalues A n are proved to be simple. Finally, efficient numerical methods 
of estimation of eigenvalues A„ are applied to obtain lower and upper numerical bounds for 
the first few eigenvalues up to 9th decimal point. 



1. Introduction 

Let (Xt), t > 0, be the one-dimensional Cauchy process, that is a one- dimensional symmetric 
a-stable process for a = 1. Let us consider the Cauchy process killed upon first exit time 
from D for D = (0, oo) and D = (—1, 1). The purpose of this article is to study the spectral 
properties of the transition semigroup of this killed process, defined by 

P t D f{x) = E x (f(X t ) ; X s G D for all s G [0, t]) , / G D>{D), 



and its infinitesimal generator Ad-, which is the operator —\J—-£? with a Dirichlet exterior 
condition (on D c ); see the Preliminaries section for a formal introduction. The key problem 
in our paper is the description of eigenfunctions and eigenvalues of Ad and P t D . The study of 
the spectral theoretic properties of the semigroups of killed symmetric a-stable processes has 
been the subject of many papers in recent years, see for example [H EJ [3j [HJ EEl EEl fT8l [19] . 
Our paper is a continuation of the work of Banuelos and Kulczycki pQ . 

In the first part of the paper (Sections [3H7|), the identification of the spectral problem for 
P t D and the so-called mixed Steklov problem in two dimensions, a method developed in pQ, 
is applied for the case of the half-line D = (0, oo). Instead of searching for a function / 
satisfying Pff(x) = e~ xt f(x) for x G D, f(x) = for x G D c , we solve the equivalent mixed 
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Steklov problem 



(1.1) 



Au(x,y) 







x e R, y > 0, 



(1.2) 
(1.3) 



d_ 

dy 



u(x, 0) 
u(x, 0) 



Au(x, 0), 



x e D 



0. 



x <£ D, 



where A = ^2+^2 is the Laplace operator in R 2 . The relation between / and u is here given 
by u(x,y) = E x f(X y ). In this way a nonlocal spectral problem for the pseudo-differential 
operator on R (or its semigroup {P t D ) on a domain D) is transformed into a local one for 
a harmonic function of two variables, with spectral parameter in the boundary conditions. 
From the point of view of stochastic processes, this corresponds to the identification of the 
jump-type process (X t ) with the trace left on the horizontal axis by the two-dimensional 
Brownian motion. Similar or related methods were also applied e.g. by DeBlassie and 
Mendez- Hernandez [TTl [THl [T9], [33] , and the idea can be traced back to the work of Spitzer [42J, 
see also [31]. 

When D = (0, 00), the spectrum of Ad is equal to (—00, 0] and is of continuous type, so 
there are no eigenfunctions of Ad in L 2 (D) (this follows easily from scaling properties of Ad] 
see also Theorem [3] below) . It turns out, however, that for all A > there exist continuous 
generalized eigenfunctions ip\ G L°°(D). More precisely, we have P^tpx = e~ xt ip\. Using the 
identification described in the previous paragraph, an explicit formula for ip\ is derived in 
Section El see fl3~T9|) and (l3~20l) . 

Surprisingly, to our knowledge, there are no earlier works concerning the spectral problem 
Pff(x) = e~ xt f(x) for x 6 D, f(x) = for x £ D c for the half-line D = (0, 00), or the 
equivalent problem (ll.ll) - (ll.3p . However, there is an extensive literature concerning the 
related sloshing problem in the half-plane, i.e. the problem given by (II. ip . (11.21) and the 
Neumann condition 



in place of the Dirichlet one ( 11.31) . The sloshing problem is one of the fundamental problems 
in the theory of linear water waves, see e.g. [23] for a historical survey. The explicit solution 
of the sloshing problem in the half-plane for D = (0, 00) was first obtained by Friedrichs and 
Lewy in 1947 [22], see also [121 [261 [29] . Both methods and results of the Section [3] are closely 
related to their counterparts for the sloshing problem in the half-plane. 

Sections [Hand [5] are rather technical and the remainder of the article relies on their results. 
Certain holomorphic functions play an important role in the derivation of and one of 
these functions is studied in Section [H In particular, the Fourier-Laplace transform of ip\ is 
derived, see ( I4.7|) . The formula for tp\ is of the form tp\(x) = sin(Ax + f ) — r(\x), where r 
is the Laplace transform of a positive integrable function. In Section [5] we obtain estimates 
of the function r. 

In Section [6] it is proved that that ip\ yield a generalized eigenfunction expansion of Ad for 
D = (0, 00) in the sense of |24J, see e.g. [35l |4"T] and the references therein for similar results 
for differential operators. In other words, the transformation Uf = (f, ipx) is an isometric (up 
to a constant) mapping of L 2 (D) onto L 2 (0, 00) which diagonalizes Ad, RAd/ = —^Ad/, 
see Theorem [31 

The spectral decomposition and results of Section H] enable us to derive an explicit formula 
for the kernel function pf(x, y) of P t D , i.e. the transition density of the Cauchy process killed 



u(x,0) = 0, 



x £ D, 
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on exiting D = (0, oo) (or the heat kernel for — y — with Dirichlet exterior condition on 

D c ), see Theorem H] in Section [TJ This extends the results of [El El EES], where two-sided 
estimates for p®(x,y) were obtained (in a more general setting). As a corollary, we obtain 
an explicit formula for the density of the distribution of the first exit time from (0, oo), see 
Theorem [5j This gives even a new result for 2-dimensional Brownian motion, see Corollary [2J 
Namely, we obtain the distribution of some local time of 2-dimensional Brownian motion 
killed at some entrance time. 

The spectral problem for the interval D = (—1, 1) is studied in the second part of the 
article (Sections IS HTTi) . We remark that due to translation invariance and scaling property 
of (Xt), the results for (—1,1) extend easily to any open interval. It is well-known that 
there is an infinite sequence of continuous eigenfunction ip n e D such that Ao<Pn = — A n </> n 
on D, ip n = on D c , where < Ai < A2 < A3 < ... — > 00. Each (p n is either symmetric 
or antisymmetric. The study of the properties of ip n and A n , dates back to the paper of 
Blumenthal and Getoor [5], where the Weyl-type asymptotic law was proved for a class of 
Markov processes in domains. In [5] (formula (3.6)) it was proved that X n /n — > 7r/2 as 
n —>■ 00. Over the last few years, there have been an increasing amount of research related 
to this topic see e.g. [H [21 ESI ESI ESI EH ED E2] and the references therein. In [TJ it was 
shown that A n < ^f. The best known estimates for general A n , namely ^f- < X n < ^ , 
were proved in [15], Example 5.1, where subordinate Brownian motions in bounded domains 
are studied. The simplicity of eigenvalues was studied in Section 5 of p], where A 2 and A 3 
are proved to be simple (simplicity of Ai is standard), and in [32], where all eigenvalues are 
proved to have at most double multiplicity. All these results are improved below. 

In Section [8] approximations (p n to eigenfunctions (p n are constructed by interpolating the 
translated eigenfunction for the half-line + x) and — x) with A = ^ — |. It is then 
shown that Ad<Pu. is nearly equal to —\(p n - This is used in Section [9] to prove that 



and that the eigenvalues A n are simple, see Theorem |6j Finally, various properties of ip n are 
shown in Section ESI see Corollaries [3J-[5j 

An application of numerical methods for estimation of eigenvalues to our problem is de- 
scribed in the last section. To get the upper bounds we use the Rayleigh-Ritz method for the 
Green operator, and for the lower bounds the Weinstein-Aronszajn method of intermediate 
problems is applied for fll.ip - fll.3l) . The numerical bounds of ca. 10-digit accuracy are given 
by formula flll.ip . 

Although probabilistic interpretation is the primary source of motivation, we use purely 
analytic arguments. In fact, the Cauchy process and related probabilistic notions are only 
used in Section[2Jto give a concise definition of the killed semigroup (P t D ), and in Appendix A. 



We begin with a brief introduction to the Cauchy process (X t ) and its relation to the Steklov 
problem. We only collect the properties used in the sequel; for a more detailed exposition the 
reader is referred to pQ or [T^l 130] . For an introduction to more general Markov processes, see 
e.g. [H [2Tj HO]- In the final part of this section, basic facts concerning the Fourier transform, 
the Hilbert transform and Paley- Wiener theorems are recalled. 




n > 1, 



2. Notation and preliminaries 
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The one-dimensional Cauchy process (X t ) is the symmetric 1-stable process, that is, the 
Levy process with one-dimensional distributions 

P x (X t G dy) = p t (y - x)dy = - * dy. 

n t 2 + (y — x) 2 

Here P x corresponds to the process starting at x G R; E a . is the expectation with respect 
to P x . Clearly, the P^-distributions of (X t + a) and (bX t ) are equal to P^+a-distribution 
of (Xt) and P ^.-distribution of (Xu) respectively; these are the translation invariance and 
scaling property mentioned in the Introduction. The transition semigroup of (Xt) is defined 
by 

P t f(x) = EJ(X t ) = f * Pt (x), f G L p (R), p G [1, oo], t > 0, 

and Pof(x) = f(x). This is a contraction semigroup on each L P (R), p G [l,oo], strongly 
continuous if p G [l,oo), and when / is continuous and bounded, then P t f converges to / 
locally uniformly as t \ 0. The infinitesimal generator A of (P t ) acting on £ 2 (R) is the 
square root of the second derivative operator. More precisely, for a smooth function / with 
compact support we have 



„ , ~f 1 [°° f(y) - f(x) 



where the integral is the Cauchy principal value. 

Throughout this article, D always denotes the interval (—1, 1) or the half-line (0, oo). The 
time of the first exit from D is defined by 

r D = inf {t > : X t i D} , 

and the semigroup of the process (Xt) killed at time td is given by 

P t D f(x) = E x (f(X t ) ; X s G D for all s G [0, t}) = E x (f(X t ) ; t < r D ) , 

where t > 0. This is again a well-defined contraction semigroup on every L P (D) space, 
p G [1, oo], strongly continuous if p G [1, oo). If / continuous and bounded in R and vanishes 
in (— oo, 0], then Pff converges to / locally uniformly as t \ 0. The semigroup (P t D ) admits 
a jointly continuous kernel function p®(x, y) (t > 0, x, y G D); clearly, p^(x,y) < p t (y — x) < 
-=j. By Ad we denote the infinitesimal generator of (P t D ) acting on L 2 (D). Since this is a 
Friedrichs extention on L 2 (D) of A restricted to the class of smooth functions supported 
in a compact subset of D, we sometimes say that Ad is the square root of Laplacian with 
Dirichlet exterior conditions (on D c ). 

Let us describe in more details the connection between the spectral problem for the semi- 
group (P®) and the mixed Steklov problem (II .ID — (II .3D . established in pQ. The main idea is 
to consider the harmonic extension u(x,y) of a function / to the upper half-plane x G R, 
y > 0. Let / G L P (R) for some p G [1, oo], and define 



u(x, y) = P y f(x) = - J ^ y2 + (z _ x)2 f(z)dz. 

Then u is harmonic in the upper half-plane R x (0, oo), and if p G [1, oo), then u(-,y) 
converges to / in L P (R) as y \ 0. Conversely, for p G (1, oo), if u(x, y) is harmonic in the 
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upper half-plane and the L P (R) norms of u(-, y) are bounded for y > 0, then u(-, y) converges 
in L P (R) to some / when y \ 0, and u(x,y) = P y f(x). By the definition, 



w(x, 0) = lim 



<9y y\o y 

pointwise for all x G R. When / is in the domain of A, then the above limit exists in £ 2 (R) 
and it is equal to Af. 

Our motivation to study the mixed Steklov problem (jl.ip — (11.3p comes from the following 
simple extension of Theorem 1.1 in [1] to the case of unbounded domains. A partial converse 
is given in the proof of Theorem [2] in Section [31 

Proposition 1. Let D = (0, oo) and A > 0. Suppose that f : R — > R is continuous and 
bounded, f(x) = for x < 0, and u(x,y) = P y f(x). If Pff(x) = e~ xt f(x) for all x G D, 
t > 0, then u satisfies (11. II) - (11.31) . 

Proof. Formulas (II .ip and (II. 3p hold true by the definition of it. Since P y fix) = e~ Xy f(x), 
we have 

u{x, y) - u{x, 0) P y f(x) - fix) e-^ - 1 ^ P y f(x) - P y D f(x) 



y y y y 

As y \ 0, the first summand converges to —Xf(x). The second one is estimated using 
formula (IA.2P from Appendix A (see also formula (2.9) in [1]). If < y < x, we have 



Pyf(x) - P y D f{x 



y 



< / — ' yK ' ' \f{z)\dz 



P y (Z-X)-p y (X,Z) 

y 

/ 1 y y \ , 2/(2 + log 



9 ' 9 ' 9 / 9 

71 n \X Z X A Z XZ A / 7TX^ 



and this tends to as y \ 0. Therefore, (11.21) is also satisfied. □ 

Finally, we briefly recall some standard facts and definitions. The Fourier transform 
of a (complex- valued) function / G £ X (R) is given by f(x) = J f(t)e~ ttx dt; this can be 
continuously extended to L P (R) whenever 1 < p < oo. For p G (1, oo), the Hilbert transform 
of / G L P (R), denoted Hf, satisfies (HfY(t) = — fit) (i sign t). This is a bounded linear 
operator on L P (R), and for almost all t, 

(2.1) Hfit) = - W I" P&ds. 

If / is Holder continuous, then the above formula holds for all t G R and Hf is continuous, 
see e.g. [43] . 

Let C + = {z G C : Imz > 0} and C + = {z G C : Imz > 0}; C_ and C are defined in 
a similar manner. Let 1 < p < oo. If F is in the (complex) Hardy space _£P(C + ), i.e. F is 
holomorphic in C + and the L P (R) norms of F(- + is) are bounded in e > 0, then, as e \ 0, 
F(- + is) converges in L P (R) to some / G L P (R), which is said to be the boundary limit of 
F . In this case 

(2.2) Im/ = #(Re/) and Re/ = -H(Imf). 
We also have 

(2.3) Hf(t) = -Hf(-t), where f{t) = f(-t). 
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The following version of Paley- Wiener theorem is important in the sequel, see e.g. [20]. For 
p G (1, oo), a function / G L P (R) is a boundary limit of some function F G H P (C + ) if and 
only if / vanishes in (— oo,0). In this case 

(2.4) F(z) = — [ f(x)e lzx dx, z G C+. 



2tt . . 

3. Spectral problem in the half-line 

Notation. To facilitate reading, in this section we strive to use the following convention. We 
use small letters to denote functions of the real variable and capital letters for functions on 
the upper half-plane C+. Real- valued functions are denoted by Greek letters, whereas Latin 
letters are used for complex-valued functions. 

We study the eigenproblem (jl.ip - fll.3B for the half-line D = (0, oo) using methods which 
were earlier applied to the sloshing problem with semi- infinite dock, see [121 1221 126] . The 
solution u is given as the imaginary part of a holomorphic function F of a complex variable 
z = x + iy, x G R, y > 0. Such a function is automatically harmonic, hence (11.1 1) is satis- 
fied. Using the Cauchy-Riemann equations, we may restate (11.21) and (jl.3p in the following 
equivalent form: 

(3.1) lm{iF'(x) + XF{x)) = x > 0, 

(3.2) ImF(x)=0 x<0. 
Observe that for all $ G R and t < 0, the bounded holomorphic functions 

F(z) = e iAz+M? , 

p(g\ ^tXz—i arctan t 

satisfy (13. ip . and for for alH > the bounded holomorphic function 

F(z) = e tXz 

satisfies (13.21) . This suggests searching a solution of the form: 

r0 

iAz+ii? / / j.\ tXz—i arctan t 



(3.3) F(z) = e lAZ+lV - / g(t)e tAz ~ t&Ictant dt, Rez > 0, Imz > 0, 

J — oo 

POO 

(3.4) F(z)= / g(t)e tXz dt, Re^ < 0, lmz>0, 



where q is an unknown real function, say in some L P (R), p G (1, oo), and d G R. The values 
of F given by (13.31) and (13. 4p must agree when Rez = 0, Imz > 0, that is, 

e ix(t) g(t)e itXy dt = e~ Xy+il3 , y > 0, 

where xif) — arctant_ = arctan(max(— t, 0)). Replacing Xy by — s yields that 

/oo 
e ix{t) g(t)e- its dt = e s+i *, s < 0. 

■oo 

The right-hand side is the Fourier transform of g{t) = t^j^- Therefore, formula (13.51) is 
equivalent to the condition: 

(3.6) the function a(t) = e* x(t) g(t) - g(t) satisfies a(s) = for s < 0. 
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Note that both g and g are in L P (R), so that a is well-defined and a G L P (R). The foregoing 
remarks can be summarized as follows: any real function g G L p (R) satisfying (13.61) yields a 
solution to the problem (I3.1I) - (I3.2I) . 

By Paley- Wiener theorem, (13.61) is satisfied if and only if a is the boundary limit of a unique 
function A in the Hardy space H P (C + ) in the upper half-plane C + = {zGC : \mz > 0}. 
Such a function A can be derived as follows. Later in this section (formula (I3.12p ; see also 
Appendix B), a function B holomorphic in C + and continuous on C + is defined, such that 
ix(t) — B(t) G R for all t G R. The function 

e- B ^a{t) = e^)-m Q{t) _ e -B{t) g{t) 

is therefore the boundary limit of e~ B ^A(z). Note that e lx ^~ B ^ is real. The function 
g{t) = T^Yfjj. is the boundary limit of a meromorphic function G(z) = ^rj^- The function 
G has a simple pole at i, so that G(z)(e~ B ^ — e~ B ^) is holomorphic in C + . It follows that 

(3.7) e- B{t) a(t) + g(t)(e~ m - e' 13 ^) = e**®-*® Q (t) - e - B ^g(t) 
is a boundary limit of 

A(z) = e~ Biz) A(z) + G(z)(e- B ^ z) - e~ B{i) ), z G C+. 

Since G and A are in H P (C+), and \e~ B ^\ is bounded (see fjM ). we must have i G H P (C+). 
Let G(z) = e 2 "^ 1 ^. z . Note that by (EZD, the boundary limit of the function A{z) -e~ B{l) G{z) 
(belonging to H P (C + )) is equal to 

(3.8) e ix(t)-B(t) _ e -B®g( t ) _ e - B ®g{t), 

which is real for all t G R. The real part of the boundary limit of an H P (C + ) function is 
the negative of the Hilbert transform of its imaginary part. Therefore, the function defined 
by (13. 8ft is the Hilbert transform of the constant 0, and so it is identically 0. It follows that 

(3.9) e *W-B(t) e (t) = e - B{l) g(t) + e - B ^g(t) = 2 Re (e~ BW (?(t)) , t G R. 

Also, A(z) — e^ B ^G(z) has a boundary limit 0, so it is identically zero in C+. Hence, for 
z G C+, 

A(z) = e B{z) (A(z) - G(z)(e- B ^ - e~ B{i) )^j 

£ -iti e B(z)-B(i) e iti 1 _ e B(z)-B(i) 

2tt 1 — iz 2tt 1 + iz 

Since |e B ^^| is bounded by a constant multiple of 1 + y\z\ (see ( 1B.8|) ) . A defined by the 
above formula is in H P (C + ) for any p G (2, oo), and g given by (I3.9P is in L P (R). Also, the 
boundary limit of A is the function a defined in ( 13.61) (this can be verified e.g. by a direct 
calculation), so that g indeed is a solution to (I3.6p . 

We now come to the construction of the function B. We want it to be holomorphic in C + 
and continuous in C + , and ix(t) — B{t) is to be real for all t G R. Therefore, 

(3.10) lmB(t) = x(t) = arctan(t_), t G R. 

Clearly B is not in H P (C + ), so that ReB(t) cannot be expressed directly as the Hilbert 
transform of lmB(t) = x(t). We can, however, apply the Hilbert transform to 1m B'(t) = 
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x'{t), which is an L 2 (R) function. It follows that 

1 f° 1 

Re B'{t) = -H{Im B') (t) = - pv / — - ds, t e R, 

n J -oo \t ~ S )\l- + S ) 

the integral on the right-hand side being the Cauchy principal value for t < 0. This equation 
is studied in Appendix B. It follows that up to an additive constant, which we choose to be 
zero, we have ReB(t) = r](t), where 77 is given by (IB. 1[) . By (1B.2P and (IB.6I) . for all t € R, 



1 /"* log \s\ 



(3.11) 



1 r i< 

B{t) = ix(t) + r]{t) = i arctan(t_) + log \/l + t 2 / - 

^ Jo 1 

, , If log It — si , 
= 1 arctan(t ) H — / — as 



+ s 2 



ds 



This formula is easily extended to complex arguments, whenever Imz > 0, we have 

1 r° log(z- s) 

7T /___ 1 + S 2 



(3.12) B{z) = -\ ° V 2 J ds, 



provided that the continuous branch of log is chosen on the upper half-plane C + (i.e. the 
principal branch with logs = log \ s\ + m for s < 0). We emphasize that (13.111) and (13.121) 
agree for z = t < (see also Section @] and Appendix B). 

For the explicit formula for g, B{i) needs to be computed. By (IB. 101) and (IB. 111) . 

ip '°f-;) rfs = lf'°f+;) ds 

J-oo 1 + S 2 7T J Q 1 + S 2 



1 f°° log(l + s 2 ) , i f°° § -arctans , log 2 ztt 
+ - / ^— cis = + — 



(3.13) 

" 27T 7 1 + S 2 ' 7T J Q 1 + S 

Now ([3i| yields that 

/V(#-f) 1 

g(t) = 2e B «-^ Re (e^Vt)) = 2e"W Re ( ^~ 

_ V2 v{t) cos(^-f)+tsin(^-f) 
2tt 1 + t 2 

Since $ G R is arbitrary, we conclude that there are two linearly independent solutions for 
g, corresponding to 1? = ~ and "d = ^ respectively, 

The solution to (I3.ip - (l3.2p corresponding to $ = | and g as above is therefore given by 

(3.14) Fiz) = e iXz+i % - — / e v(t) t\z-i octant df R > j > 

1 ^ 1 ^ 2n J^l + t 2 - ' 

(3.15) F(z) = — -e^e^dt, Rez<0, lmz>0. 

By (IB. 51) . we have g G L 1 (R), and so F is bounded and continuous. Furthermore, it can 
be easily verified that the solution corresponding to $ = ^ and g is given by — M. Since 
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10 12 



(b) 



Figure 1. (a) Graph of ipx\ (b) Graph of the remainder term r(x) = sin(x + 

i)-V>iO) 

decays at infinity as \t\~^ 7 F'(z) has a singularity of order \z\~^ at zero and it is not bounded 
near 0. For that reason, in the sequel we only study the solution F(z) given by (I3.14p - (l3.15p . 
Since e - iarctan * = (l-it)/V±T~P and = e^ - *) VTTt 2 (see (El])), we can rewrite (JXHD 

as 

V2 f°° 1 + it 
2% 



(3.16) 



F(z) 



- V (t) e -tXz dt 



Rez > 0, Imz > 0. 



1 + t 2 

Therefore, we proved the following theorem. 

Theorem 1. The bounded solution of (II. II) - (11.31) for D = (0, oo) is given by 

u(x, y) = e~ Xy sin(Aa: + |) 
(3-17) ,/2 t f°° tcos(tXy) — sin(tXy) /If 



271 



for x > 0, y > 0, and 
(3.18) 
for x < and y > . 



u(x,y) = ^ 



1+t 2 

sin(t\y) 



log(t + s) 



1 + s 2 



ds 1 e~ tXx dt 



i + t 2 exp u 



logg + g) 
1 + s 2 



rfs ) e tXx dt 



The main result of this section, stated below, follows from Theorem[T]and a partial converse 
to Proposition [TJ 



Theorem 2. Let D = (0, oo). For A > 0, the function 
(3.19) if) X (x) = sin(Ax + f)-r A (x), 

w /i ere 

V2 r t (I f°°log(t + s 

to 1 — 6XP '^ 1 



(3.20) 



r A (x) 



r(Ax) 



1 + t 2 



7T 



1 + ^ 



x > 0, 



-ds e 



-tA:c 



dt, 



is the eigenfunction of the semigroup (P t ) acting on C(D) corresponding to eigenvalue A. 

Proof. With the notation of Theorem [TJ we have ip\(x) = u(x,0); we extend ip\ to be on 
(— oo, 0]. Since u is harmonic and bounded in the upper half-plane, we have P y ip\(x) = u(x, y) 
(y > 0, x G R). Since u satisfies (11.21) . for all x > 0, -(P y il)\(x) — ip\(x)) converges to 
— Xtp\(x) as y \ 0. We will now prove (formula (13.261) ) that this convergence is dominated 
by an appropriate function. 
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Below we assume that A > 0, x > and < y < v. By formula (13.171) . we have 
(3.21) V 



\/2 f°° tcos(t\y) - sin(tAy) - e~ Xy t 



2tt 



- V (t) e -t\x dt 



Since |1 — cosz\ < \z — sinz\ < y, |1 — 2 — e 2 | < ^- and Xy < 1, we have 

|f cosCfAy) - sin(tAy) - e- A ^| < A 2 t ^ + ^ + y 2 < A 2 t(l + t 2 )y 2 . 
Using also < e c ^(l + t 2 )"i (see flB~4|) ). and then (1 + 1 2 ) > t 2 , we obtain 



(3.22) 

Furthermore 
(3.23) 



i tcos(t\y) — sm(t\y) — e Xy t 



(l + t 2 )y 



- v (t) e -tXx dt 



< e^\ 2 y 



t(l + r 



(1 + t 2 ) 



e- tXx dt < e*\ 2 y / Vte~ tXx dt. 



i_ 

Ay 



Vte~ tXx dt < (Ay)-t / H e _tAx dt < 



r( 



(X 2 xy) * 

In a similar manner, but using |t cos(t\y) — sin(tAy) — e^ Ay t | < t + tAy + t < 3t, we obtain 
tcaa(t\y) - sin(tAy) - e~ Aj/ t 



(3.24) 

and 
(3.25) 



[l + t 2 )y 



- V (t) e -t\x dt 



< 3e» 



t , 3e' 

-^e~ tXx dt < — 



r - e -tXx dt 



t—*e- tXx dt < (\y)i I r^e~ txx dt 



4„-a* r(i)(Ay)f 



(Ax 



Formulas (j3.2ip - (13.25p yield, after simplification, that 



(3.26) 



< 



2v / 2eiT(|)A5 2/ i Cl (\) y i 

3 

7TX4 



3 

;r 4 



with some constant ci(A). 

We are now going to replace P y by P® in (13.261) . In Section [5] it is proved (using only the 
definition (EDS]) and (I3T2T?]) of ip x ) that |^a(»| = |^i(Ax)| < 2V'Aa;, see (Oil . This and ([Oil 
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yield that for < y < x we have 



11 



y 



< 



• Py(Z ~ X) -Py(Z~ X) 



y 



\if>x(z)\dz 



(3.27) 



< 



2y/\ 



71 



miu ( — , ) y/zdz 



When < x < y, in a similar manner 



8y / A|/(3 x /i- jy) 8V\y 8y/\y* 

o 2 ° 

6ixx z 



3 

7TX2 



3 

7TX4 



(3.28) 



< 



2y/X 



7T 



16 16^2/3 



Q 2 
37TX4 



Finally, by ( 13. 26ft . ( 13.271) and ( 13.28j) . there is a constant 02(A) such that 



(3.29) 



< 



c 2 (A)y 
3 

;r 4 



For any fixed x > and i > 0, the one-sided derivative of e xt P t D ip\(x) with respect to t 



equals 



_d_ 

dtl 



( xt pD I / \\ , e^P t D +y M*)-e xt P t D Mx) 

(e P. ip\(x)) = hm 

y\o y 

= lW<-> fp°(,,,) P > (z) - e " >,fe(z) 
y\° Jo y 



dz. 



Since p?(x, z) <p t (z — x) < |, by ( 13. 29ft we have 



P* (x, z) 



f/2 



< 



1 — Pt (x,^)dz 



24 



<c 2 (A)^^ p?{x,z)dz+^^ < (l + ^c 2 (A)yi 

The right hand side tends to zero as y \ 0, so that 



for allt > and x > 0. By (13.291) (with both sides multiplied by e Xy ), this also holds for t = 0. 
Finally, the function e xt P t D ip\(x) is continuous with respect to t for each x > (this follows 
from weak continuity of pf{x,z)dz with respect to t, which is a consequence of stochastic 
continuity of the killed Cauchy process; one can also prove this using the explicit formula for 
p t and ( 1A.2I) ). It follows that e xt P t D tl>\(x) is constant in t > 0, and since Pqi})\(x) = ip\(x), 
this completes the proof. □ 
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Remark 1. Since r(x) > 0, the functions ip\ are clearly not in L 2 (D), so the above result 
does not provide any information about the L 2 (D) properties of the operators Pf . This 
problem is studied in Section [6j 

Remark 2. Also the derivatives ip' x are locally integrable eigenfunctions of P t , continuous 
in D, but not at 0. As this is not used in the sequel, we omit the proof. 

Remark 3. The functions tp\ can be effectively computed by numerical integration. In- 
deed, (jR2]h flR6l) and the identity 

/"* log U , , % / t . / . \ T • / ■ w 

/ -du = arctan t log t H — (Li 2 Ut) — L12I— it)), 

Jo 1 + u 2 2 

where Li 2 (z) is the dilogarithm function, yield that 

^ A (x) = sin(Ax + |) - — / ^exp — (Li 2 (it) - Li 2 (-i*)) e' tXx dt. 

^ JO (1 +t 2 )4 V 271 " / 

4. Properties of the function B 

In this section we study the properties of the function B. As an interesting corollary, the 
Laplace transform of the eigenfunctions ip\ is computed. 

The function B defined by (I3.12p extends to a holomorphic function on C \ (— 00, 0], 
satisfying B(z) = B(z) (see also Appendix B). Therefore B is defined on whole C, it is 
holomorphic in C \ (—00, 0] with a branch cut on (—00, 0], and it is continuous on C + . The 
following properties of B will play an important role. 

When Imz > 0, we have 

/ s 1 f° tog(z-s) , 1 f° log(-z-s) , 
B(z)+B(-z) = - r. o ds + - S : , ' ds 



J -00 1 + S 2 7T J_ OQ 1 + S 2 

1 rio& +S ) d 1 f M i± 4 zi£(is 



IT J Q 1 + S 2 7T ./ N I + -> 

1 f°° \og(z - s) , m 
' ds . 



nj^ 1 + s 2 2 

On the right-hand side, the function s i— > log(2; — s), holomorphic (and therefore harmonic) 
in C_, is integrated against the Poisson kernel of the lower half-plane pi(s) = ^y+s 5 " The 
result is the value of log(2; — s) at s = —i. It follows that 

ITT 

B(z) + B(-z) = \og(z - (-i)) - - = log(l - iz). 



By B{z) = B(z) we get B(z) + B{—z) = log(l + iz) whenever Imz < 0, and so 
(4.1) e B{z) = (l-iza(z))e- B{ - z \ 

where a(z) = 1 when I1112 > and a(z) = —1 when Imz < 0. A similar relation for 77 was 
used earlier in (13.161) . see also ( IB. 31) . By continuity of B(z) in C + the formula (14.11) is also 
valid for z G R if we let a(z) = 1 for z < and a(z) = —1 for z > 0. For completeness, we 
let ct(0) = 0. 
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By f )3.14p . (13.201) . the relation between F, ip\ and r A , and using B(t) = r](t) + iarctant_, 

a/2 f°° e 13 ^ 

(4.2) r(x) = j— r(t)e~ tx dt, where r(t) = Im -. 

2ir J 1 + V 

Note that by the definition, r(t) = for t < 0. In the sequel, we need the Hilbert transform 
of r, which can be computed as follows. The function e B ^ /(l + z 2 ) is meromorphic in the 
upper half-plane with a simple pole at i, so that the function 

e B(z) 1 e B(i) x e B(ij 

1 + z 2 ~~ 21 + iz ~ 2 1 -iz 
is holomorphic in C+. In fact it is in H P (C + ) for p G (1, oo), see (IB. 81) . Its boundary limit 
on R is equal to 

e B(t) 1 e B(i) 1 e Wj e B(t) _ y2cos f - ty/2 sin f 



1 + t 2 21 + zt 21 -it l + t 2 

and the imaginary part of this function is just t(— £). Therefore, the Hilbert transform of 
r(— t) is the negative of the real part of the above function. It follows by (12.31) that for t G R, 

. . /N e B ^ - v^cosf +tv^sinf 
(4.3) Hr(t) = Re ^ §- . 

We are now able to compute the Laplace transform Cipx of By a direct computation, 
we have 



f°° i cos f + t sin f 
(4.4) y sin(x + f )e~ tx dx = — §■ , t > 0. 

On the other hand, by Fubini's theorem and (14.21) . 



oc 







r(x)e- tx dx = — ri^-ds = --Hr(-t), t > 0. 

V ; 2tt J t + s 2 y h 

By (14. 3 p we have 

, s f°° , ^ tn. V2 e B ^ cosf + tsinf 
4.5 / r (x )e~ tx dx = + ? — , t > 0. 

In particular, 

(4.6) / r(x)dx = cos . 

1 ^ io 8 2 

Formulas (14.41) and H4.5|) give 

/"°° \/2 e B W 

CMt) = J Mx)e~ tx dx = ^YT¥ ' t>0 ' 

By scaling and the uniqueness of the holomorphic continuation, we obtain the following 
result. 

Corollary 1. The Laplace transform ofip\ is equal to 

f°° \/2 \e B ^ 

(4.7) Wx{z)= i, x {x)e- zx dx=^-—— I) Re^>0, 

Jo 2 y + z 2 

where B(z) is given by (I3.12p . 
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5. The remainder term 

This section is devoted to a detailed analysis of the remainder term r\, see (13.201) . Recall 
that r\(x) = r(Xx), where 

Since r is the Laplace transform of a positive function, it is totally monotone, i.e. all functions 
(— iyM n ) are nonnegative and monotonically decreasing. In most of our estimates we simply 
use the inequality — C < jf^ds < for t > and formula (IB. 12ft from Appendix B. The 
L 1 (R) norm of r, however, we already calculated in (14.61) . 

/*°° 7T /2 

(5.2) / r(x)dx = cos — — —— G (0.216, 0.217). 

Jo 8 2 

Since -|- < 7-7^, by Fubini's theorem, 

t+s — 2vts' J ' 

, , „ 3l 1 f°° f°° t S 1 , , 

(r{x)j ax < — - / / - g- aids 



2vr 2 J J (l + tz)! (l + s i)zt + s~ 



<5 ' 3) <JLfr^*v=-a<o.o3, 

-4vr2W (l + t 2)f J 7T(r(i)) 2 



In a similar manner, t— > , 7= a , so that 



2C / .„ \ 2 2C 



(5.4) / (r{ x )Ydx> e —[ f t t - dt) =^f> 0.012. 



2tH (1+t*)* y 9tt 

For a; > 0, we have 



(5.5) 



V2 r . _ to . v/2 



< — / te~ tx dt < 



2ir In 2tyx 2 



In a similar manner, 



(5.6) (-l) B rW(x)< V2!/ ' 1 l: 



2tt s"+ 2 



Also, 

(5.7) r(x) < r(0) = sin f = ^ 2 ~ ^ < 0.383, 

and 



CAUCHY PROCESS ON HALF-LINE AND INTERVAL 15 

It follows that for x > 0, 

l^i 0*0 1 ^ l sin ( x + f ) - sin 1 1 + \ r ( x ) - r (°)l 

<x+ / r (y) lay < x + \ — . 

J0 V 7T 

Since clearly |'0i(^)| < I sin(x + f )| + \r(x)\ < 2, we have 

(5.9) < min hr + y^> 2 j < min (2v^, 2) . 

This property was already used in Section [3] in the proof of Theorem [2j 
Finally, since r' < 0, the first zero of ip' x is greater than It follows that 



y/2 f°° te 



(5.10) ||^<l+r A (£) = l + r(§)<l + ^- / - T dt<1.14 



2tt 



o (1 + t 2 ) 



4 



for the last inequality, integrate by parts the left hand side of formula 3.387(7) in [25]. The 
estimate (15.101) is only used in Corollary O where a weaker version of (I5.10P would only 
result in a larger constant in (110.21) . In fact, for the present constant 3, we only need that 
IIV'aII < 1.19, which is easily obtained by estimating (1 + 1 2 )~ 5//4 in the integrand in (15.101) 
by a constant on each of the intervals [|, ^p] with k — 0, 1, 7, and [4, oo). 

6. Spectral representation of the transition semigroup for the half-line 

Let D = (0, oo). In this section we study the L 2 (D) properties of the operators Pf . For 
/ G C C (D), define 

POO 

(6.1) n/(x)= / f(\)ij)\(x)d\, xeD, 

Jo 

where ip\ = sin(Ax + | ) — r\(x) is given by (13 . 19[) and (13.201) . Note that 

POO 

F 1 (x)= /(A) sin(Ax + f )d\, x E D, 

Jo 

satisfies \\Fi\\ 2 < C\ ||/|| 2 - Also, for 

POO 

F 2 (x) = / f(\)r x (x)d\, x E D, 

Jo 

we may apply (I5.5P and (15.71) to obtain 



oo poo poo 



(F 2 (x)fdx< / / / \f(fi)\\f(X)\r(fix)r(Xx)dfidXdx 
Jo Jo Jo 

poo poo / poo \ 

I U F^ 7 /Mll/Wl¥> 

poo POO / POO \ 



oo />oo 



o Jo V + X 
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which is bounded by ||/|| 2 by Hardy-Hilbert's inequality. It follows that ||n/|| 2 = 
II -Pi + -^2 1| 2 < c 3 II /II 2' anc ^ therefore II can be continuously extended to a unique bounded 
linear operator on L 2 (D). 

For / G C C (D), pP(x,y)f(\)ip\(y) is integrable in (y, A) G D x D, so that by Theorem^ 



(6.2) irn/(x) = / e- A 7(A)VA^)^A, x £ D. 

Jo 

Let f,gE C C (D) and define / fc (A) = e" fcAi /(A), # fc (A) = e" fcA *c/(A). From (Q it follows that 
PfUfk = Hfk+i and P®Hg k = ILj fc+1 . Since the operators P/ 3 are self-adjoint, we have 

poo /*oo 

/ n/(x)n^(x)rfx = / P t D Uf^(x)Ug(x)dx 
Jo Jo 

poo poo 

= / Uf^(x)P t D Ug(x)dx = / n/_ x (x)n^(x)dx. 
Jo Jo 

By induction, 

/•OO POO 

/ n/(x)Il5((x)c/x = / Uf^ k (x)Ug k (x)dx. 
Jo Jo 

Suppose that supp / C (0, A ) and supp g C (A , oo). Then we have 

POO 

Uf(x)Ug(x)dx= / U(e- kXot U)(x)U(e kXot g k )(x)dx. 



Both e~ ot f- k and e kX(>t gk tend to zero uniformly as k — > oo, and so Il(e~ fcAo '/_fc) and 
n(e fcA °*(7fc) converge to zero in L 2 (D). We conclude that 11/ and lip are orthogonal in 
L 2 (D). By an approximation argument, this is true for any f,g G L 2 (D), provided that 
/(A) = for A > A and #(A) = for A < A . 
Define 

POO 

H(A) = / {Yil A (x)) 2 dx, AC D. 

Jo 

Clearly 

/i(A)<ca||l A ||* = c 3 |A|, A CD. 

Whenever A C (0, Ao) and £> C (Ao, oo), we have 

POO POO POO 

/i(AUB)= (m A (x)) 2 dx+ (m B (x)) 2 dx + 2 m A (x)m B (x)dx 
Jo Jo Jo 

= fi(A)+fi(B). 

Finally, when A = |J^=i Ai> where C 4 C ... and |A| < oo, the sequence l An converges 
in L 2 (D) to 1 A as n — > oo. Hence m An converges to ITl^ in L 2 (D), and so fi(A) = 
Hm^oo fi(An). It follows that /i is an absolutely continuous measure on (0, oo). By an 
approximation argument, we have 

POO POO 

/ Uf(x)Ug(x)dx = / f(X)g(X)ii(dX) 
Jo Jo 

for any f,g E L 2 (D). 

Note that ip\(qx) = ipx/q( x )-> an d therefore Hf q (x) = qllf(qx), where f q (x) = /(-). It 
follows that n(qA) = qfi(A) and so /i must be a multiple of the Lebesgue measure on (0, oo), 
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say n(A) = C4\A\. This result is a version of Plancherel's theorem, where Fourier transform 
is replaced by II: 



' /»oo 

Uf{x)Ug{x)dx = c 4 / f(X)g(X)d\ 
Jo 



for any f,ge L 2 (D). 

The constant c 4 can be determined by considering /(A) = -^=l[ 1)1+ q](A), q > 0. We then 
have || /|| 2 = 1. On the other hand, 

•l+g 

r(Xx)dX. 



Uf(x) = -i- (cos(a; + f ) - cos((l + q)x + f )) - ±- I 
x^fq y/qJi 



The L 2 (D) norm of the first summand converges to as q \ 0, just as in the case of the 
Fourier sine transform. The second summand is bounded by y/qr{x) and so it converges to 
zero in L 2 (D). It follows that C4 = |. The Plancherel's theorem can be therefore written as 

(6.3) / Tlf(x)Ug(x)dx = - / f(X)g(X)d\ 



In particular, yf II is an isometry on L 2 (D). Since ip\(x) = ip x (X), for f,g G C C (D) (and 
therefore for any /, g G L 2 (D)) we also have 

n/(z)0(x)tfc = / /(A)n 5 (A)dA, 



which combined with (16. 3p yields that II 2 / = |/. We collect the above results in the 
following theorem. 

Theorem 3. The operator y/f n : L 2 (D) -> L 2 (L>) owes a spectral representation of Ad 
and the semigroup {P t D ), acting on L 2 (D), where D = (0, 00); t/jat is, /or any / G L 2 (D), 

(a) ||/|| 2 = a/| ||n/|| (Plancherel's theorem); 

(b) UP t D f(X) = e- xt Uf(X); 

(c) / is m i/ie domain of Ad if and only if XUf(X) is square integrable; 

(d) IL4 D /(A) = -AII/(A). 
Furthermore, U 2 = | Id (inversion formula). 

7. Transition density for the half-line 

The aim of this section is to compute an explicit formula of the transition density pf(x,y) 
of the Cauchy process killed on exiting a half-line D = (0, 00), or the heat kernel for Ad- 
Let us note that the transition density of the Brownian motion killed on exiting a half-line 

(0,oo) equals -J=e~~ ^ i=e~ ' 2"' , which follows from the reflection principle. For 

the Cauchy process we cannot use the reflection principle and the computation of p®(x,y) 
requires using much more complicated methods. 

Theorem 4. For D = (0, 00) and any g G LP , p G [1, 00] , we have 

POO 

(7.1) P t D g(x) = p?(x,y)g(y)dy, t,x>0, 
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where 

(7-2) V? (x, y) = V2 - - f f{ : )fi ^ f) ds, t,x,y> 0, 



and 



7.3 f(s) = -—- -J exp - / , s > 0. 



Remark 4. Note that for s > 0, / is positive continuous and bounded. This follows by the 
fact that f(s) = ^ix^ e an d f|B.4j) . The function j/) can be effectively computed by 
numerical integration. Indeed, by the same arguments as in Remark [3] we have 



Is 1 



arctan(s) 



—i 



= ~Ti ^ eXP ^-( Li 2( is ) - Li 2(-^)) 

where Li2 is the dilogarithm function. 

Proof of Theorem^ For <? G C C (D) we have IIP t D ^(A) = e~ A 'rig(A) (see ( 16.1 p and Theo- 
rem |3J). Applying II -1 = -IT to both sides of this identity yields 

P t D ^(x) = - / e- xt Ug(\)ij x (x)d\ = - / e- xt ^ x {x)^ x {y)g{y)dyd\. 
^ Jo 71 Jo Jo 

By the Fubini's theorem, (17. ip holds with 

2 f°° 

(7.4) p?(x,y) = - Mx)My)e- Xt d\. 

By an approximation argument, (17.11) holds for g G L P (R) with any p G [1, oo]. We will now 
prove (I7.2p . 

Suppose first that x < y, and let t = t 1 + t 2 , ^i,t 2 > 0. By Plancherel's theorem and 
identities ip\(x) = ip x (X), Cip y (z) = Cip y (z), we have 

p?(x,y) = - / ^ x We- tlX )mX)e- t2X )d\ 

(7.5) 

= — Ci/j x (t 1 +is)Cipy(t 2 -is)ds = — / R(z)dz, 

n J-oo ^ m Jti-ioo 

where (see (14.71) ) 

2 1 xyexp(B(i) + B(^)) 

R{z) = -Cip x (z)C^ y {t -z) = - ———JL—— — — — . 
7r ir [x z + z 2 ){y 2 + [t — z) 2 ) 

Note that R is defined on C and it is meromorphic in C\ ((— oo, 0] U [t, oo)) with simple poles 
at ±ix and t ± iy. Let z G C \ [0,t]. By (14.11) and the identity (1 + izo~(z))(l — izcr(z)) = 
1 + (zcr(z)) 2 = 1 + z 2 , we have for all z G C, 



i?(z) 



1 (1 ~ (1 - gg^g)) exp(-g(-f ) - g(-g)) 

^(1 + ^(1 + ^) 
! exp(-£(-f ) - gHg)) 

7T 0^/(1 +i|<7(f))(l+i^(^)) ' 
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Since a{^f) = for z E C \ [0,t], it follows that for z e C \ [0,t], 



t—z\\ / 7 t—z 



We therefore have R(z) = R\{z) + -fi^- 2 ) for zgC \ [0, i], where, again using (14. ip . 

exp(-S(-f)-B(-*=2)) | 



(7.6) 



and 



exp(S(f ) - B(-*=*)) £ 



2 • 



(7.7) R x {z) 



exp(S(f) + 5(S=5)) 



3/ v y 

Also, in a similar manner, 

exp(-fl(-f) 



(7.8) R 2 (z) 



y " y 



and 

(7.9) ifc(z) = 



2/ 7/ y 



The only zero of - + — is z = < 0. Hence R±(z) is holomorphic in the set {Re z > 0} \ 
[0, t] (by (I7.6P ). bounded in the neighborhood of [0, t], and it decays as \z\~ 2 at infinity 
(by (1B.9I) ). Also, R2(z) is meromorphic in the set {Re z <t}\ [0,t] (by ( 17.81) ) with a simple 
pole at bounded near [0, t], and it decays as \z\~ 2 at infinity. 
For n = 1, 2, ... let 7 be the positively oriented contour consisting of: 

• two vertical segments 71 = \t\ — ni, t± — ^1 , 75 = [ti + £1 + ml , 

• two horizontal segments 72 = [*i - -£] , 74 = [~, ti + , 

• two semi-cirles 73 = {|^| = ~, Rez < 0} and 76 = {|z — t x \ = n, Kez < ti}. 

Clearly, / ,_, R 2 (z)dz converges to J^-^ R2{z)dz as n — > 00. The integrals over 73 and 76 
converge to zero by the properties of R 2 . Finally, by (17. 9p . 

p , u T exp(B(f) + B(^)) *g /l 1 \ 

RoiZjdz — > / 7 7 ^ • 72 US 

Jo ™/(f • '7 ' • '7/ 

41 2iri/(f)/(^) 

Therefore, by the residue theorem, 

I rh+ioo r ti f(±)f(t = s\ 

(7.10) — / R 2 (z)dz = - r'L\ ds + Res(R 2 ^ y ). 

Z7ri Jti-ioo JO x y\x>~^~) 



72U74 
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In a similar manner, using (17. 7p and analogous contours 7 consisting of two segments of the 
line Re z — t±, two segments parallel to [ti,t], and two semi-circles centered at t (the small 
one) and t\ (the large one), both contained in {Re z > t\}, we obtain that 



(7.11) / R 1 (z)dz = - — j f-^-ds. 

Therefore, (f73]l . flTTTUD and flTTTTD yield that 

x ft ffs)fft=S.) 

P?(x,y) = -— / J r'\_ y s ds + Res(R 2 
xy Jo ~ + -r 



x-y 

V 



For z = we have - = — — = — — . Therefore, by 07.81) we get 

x—y x x—y y 1 J 1 ' ° 



Res(i? 2 , ^ 



1 ~ 

tx \ _ x x-y 



n(y-x) l + j^j Trt 2 + (x-y) 2 



and (17.21) follows for x < y. 

When x > y, simply note that pf(x,y) = p^(y,x) (see (17.41) or e.g. [Hj, Theorem 2.4), 
and that the right-hand side of (17.21) has the same symmetry property (this follows by a 
substitution s — t — v). Finally, for x = y simply use the continuity of p®(x, y) and /. □ 

For the next result, we need the following simple observation, similar to the derivation 
of (14.31) . By (14.11) we have Ime _B ^ s ' = — j^e B ^ for s > 0. Hence the function / defined 
by (17. 3j) satisfies 

/(*) = -T^^ = -T^-2 eB(s) = --Ime-*(-*), 8 > 0. 

71 1 + S Z 71 1 + S z 71 

If we extend / by f(s) = for s < 0, then /(— s) = - Ime" B ' s ^ for all real s. Since e~ B ^ is 
in H P (C+) for p G (2, oo) (see ( IB. 81) ). the Hilbert transform of / is given by (see (12. 31) ) 

Hf(s) = — Ree- B{ ~ s) = --e^ { ~ s) , s G R. 

71 71 

It follows that 

(7.12) #/(_ s ) = _J__^_L s>0 , and H f(0) = --. 

TV 2 1 + S 2 j(S) 71 

Theorem 5. For D = (0, oo), we have 

1 x fir 00 loef- + w) \ 

(7.13) p. (rDe(i() = __e xp y 

Using the function / defined in (17.31) . we have P x (t d G eft) = \f(^)dt. 
Proof. By Theorem H] we have 

pf(x,y)dy 

(7J n i r / , r r 1 , , 

/ / — : , t -,s dyds. 



tt Jo t 2 + (x- y) 2 J J xy § + * * 
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By a substitution w = (t — s)/y we obtain 

i /(!)/( 



y_ 



/(f) 








/(f) 




s 


Vo w 



= — — / — ^— ; — ^ 

fH 



X 



+ w 



dw 



The right-hand side equals f /(f )(-Hf(0) + #/(-§)). This, (17321) and iTH give 

1 Z" 00 t /"* f( £ ) 1 /"* x 

By substitution t> = x — y in the first integral and v = xt/ s in the third one, 

i r x t r* f( £ ) i f°° t 



s 

The result follows by differentiation and ( 17. 3ft . □ 

Remark 5. Theorem [5] can also be obtained in a more explicit manner. In fact, by scaling 
properties of X t , we have P x {td > t) = g{j ) for some function g continuous in R, vanishing 
on (—oo,0]. Furthermore, P x (td > t) satisfies the heat equation in D, i.e. ^ P x (td > 

*) = ~y Px( t d > t). For t = 1 this gives -xg'(x) = -\J g(x), x > 0. Since 



d 2 

cix 2 



(7 = Hg', we have Hg'{s) = —sg'(s) for s > 0. 



Let h(s) = -g'(— -), /i(0) = 0. Then it can be shown that /i is continuous on R, and by the 
definition flUQ) of the Hilbert transform, #7i(s) = -±##'(-1). It follows that h(s) = sHh(s) 
for s < 0, and h(s) = for s > 0. Therefore Hh — ih is a boundary limit of some holomorphic 
function in C+, and (Hh(s) — z/i(s))e~ iarctalls_ is real for all s G R. This problem can be 
solved using the method applied in Section [3j and the solution is Hh(s) — ih(s) = ce~ B( - s ^ 
with some c G R. Therefore g'(s) = —^h(—^) = - Ime - ^ 1 /^ for s > 0, and finally 

P x (r D G dt) = -f </(f) = -f Ime- B ^ x \ 

which agrees with (17.131) when c — -. The details of this alternative argument are left to 
the interested reader. 

Remark 6. The integral of p® (x, y) with respect to t G (0, oo) is the Green function of (X t ) 
on the half-line, given by the well-known explicit formula of M. Riesz, see e.g. [6]. Also, 
the distribution of X(td) (and even the joint distribution of td and X{td)) is determined 
by pf{x,y), see [27]. Explicit formulas for Green functions and exit distributions for some 
related processes in half- lines and intervals were found recently in [101 [IT] . 

Theorem [5] implies a new result for the 2-dimensional Brownian motion. Namely we 
obtain the distribution of some local time of the 2-dimensional Brownian motion killed at 
some entrance time. For the 1-dimensional Brownian motion similar results were widely 
studied and are usually called Ray-Knight theorems 
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Corollary 2. Let B t = (B^\b^) be the 2- dimensional Brownian motion and L(t) = 
lim £ _^ + y Jq X(-e,£)(b1 2 ^) ds be the local time of B t on the line (—00,00) x {0}. Let A = 
(—00, 0] x {0} and T A = inf {t > : B(t) G A} be the first entrance time for A. Then for 
any x > we have 

P^\L{T A ) G dt) = -^— 2 exp ( 1 ^ l ° g 1 ( « + a 8) 
For (x, y) G R 2 , y 7^ and t > we have 



1 /"° 

p(*>y)(L(T A ) <t) = - \ 

7T J-oc 



I v I 



y 2 + (x — m) 2 



Proof. Let r/ t = inf{s > : L(s) > t} be the inverse of the local time L(t). It is well known 
(see e.g. [42]) that the 1-dimensional Cauchy process X t can be identified with B^'^rjt). 
With this relation, we have L(T A ) = T(o,oo), where T(o,oo) = inf {t > : X t ^ (0, 00)}. This 
and Theorem give the first equality. The second equality follows by the harmonicity of 
(x,y) -> P^y\L(T A ) < t) in {(x,y) G R 2 : y > 0} and in {(x,y) G R 2 : u < 0}. □ 



8. Approximation to eigenfunctions on the interval 

In this section the interval D = (— 1, 1) is studied. Let n be a positive integer and \i n = ^ — | . 
Our goal is to show that \i n is close to A n , the n-th eigenvalue of the semigroup (Pf). 

Let q be the function equal to on (00, — ~) and to 1 on (|, 00), defined by (IC.ip in Ap- 
pendix C. We construct approximations to eigenfunctions of (P t D ) by combining the eigen- 
functions V ; Ai„(l + x ) an d Vw(l — x ) f° r half-line, studied in Section [3J For a symmetric 
eigenfunction, when n is odd, let 

<Pn(x) = g(-x)^ Mn (l + x) + g(x)^ n (l - x) 

(8.1) ^ 

= (-1) 2 cos{ii n x)l D {x) - q{-x)r^{l+x) - q{x)r tJ , n {l-x). 

For an antisymmetric eigenfunction, when n is even, we define 
<Pn{x) = q(-x)ip^(l + x) - q(x)ipn n (l - x) 

(8.2) n 

= (-1)7 sin(/i n x)l D (x) - g(-x)r /tn (l + x) + g(x)r Atn (l - x). 
Lemma 1. With the above definitions, 



1 , -1 ' 8-00 13.66 1 
(8.3) ^ n + ^n a < WI.2I+ + — . 

Proof. Note that we have 

<p n (x) - ^„(1 +x) = -(1 - g(-ar))^ Mn (l + x) - (-l) n g(x)^ n (l - x) 
= -ff(aO(Vw(l + *) + - *)) 

= <z(x)(r Mn (l + x) + (-!)%„(! - x)) - sin(/x n (l + x) + |)l[i )0 o)(x). 
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Denote h(x) = sin(/i n (l + x) + f)l[i )0O )(^) and f(x) = r Mri (l +x) + (-l) n r /in (l - x), g(x) 
g(x)/(x). By (EES), (EE!) and (ED, 

M = snp |/(o;)|<r(^)+r(^)<i^|, 



xg( 3 , 3 ) 

M 1= snp | / '(x)|<- yUn r'(^)- / i n r'(^)< 2l ' >v/2 



xe(-i,i) 



M 2 = snp ir(*)|</iy'(^)+/^%)< 



ie( 3 , 3 ) 



647T//2 ' 

4131\/2 
2567r/x^ ' 







1 

l^n Jo 



\f(x)\dx< I r fln (l + x)dx + / r lJLn {l-x)dx 
7T V2\ 1 



r(y)dy = cos • 



2 / ju „ 



notation here corresponds to that of Appendix C. By ( 1C.2I) and ( 1C.3I) . 



(8.4) 
(8.5) 



, „ . 0.605 0.156 
\A D g(z)\ < —^- + , 

4.444 0.622 



\A D g(z)\ < 



+ 



2 1 ' 



Furthermore, \g(z)\ = for z G (—1, — |) and 
(8.6) \v n g(z)\ < fM < 

Finally, for z < we have 



3 - 
0.317 



*e(-l,0). 
zg (-±0). 



7T 



sin(/x n (l + x) + |) 



fix 



x — z) 



< 



so that 
(8.7) 



7r(l — z) 1 



AN 1 ,/ \ I 0.180 
/in 

an 1 ,. \ i 0.319 



|sin(/i n (l + x) + |) I dx 



7T/X n (l - z) 



2 - 



ze (-|,o). 



(8.8) 

Since for z G (—1, 0) we have 

|./4z>£»Cz) +Hn<P n (z)\ < \(-&)*h(z)\ + |(-A)^(z)| + 
estimates ljOl- fETB|) yield that 



(8.9) 
(8.10) 



. 0.605 0.336 

\A D ip n (z) + /VyWJl < — — H 

4.444 1.258 



(-1, 



3/' 



+ 



ZG(-|,0). 
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By symmetry, estimates similar to (18.91) and f 1 8 . 1 j) hold for z £ (0, 1). The estimate (18.31) 
follows. □ 

The estimate of the L 2 (D) norm of (f n plays an important role in the sequel. We have 



:$..u) \li- — <\m\ a <Ji + — . 

V l*n 

Indeed, the lower bound follows by (15. 2p . (18. ip . (18.21) and symmetry, 



\¥n\\\ > J (sin(/i„(x + 1) + \)) 2 dx 



-4 / \q{-x)r^ n (x + l)sin(/x n (x + 1) + f) eta: 



-1 



> 1 + — cos-- — 

4/i n J ii n \ 8 2 



In a similar manner, using also (15.31) . 



V2\ 4 / 7T v 7 ^' 



r»ll2 < I 1 J I COS o - -7T I + 4 / (H//„(-M 1))) L V'- 

<| 1 + ^|Wcos*-^| + W 



4/i„y /i n ^ 8 2 y 7r(r(i)) 2 /i n " 

9. Simplicity of eigenvalues for the interval 

We continue denoting by <pj the eigenf unctions of (P/ 3 ), by Aj (Aj > 0) the corresponding 
eigenvalues, and by y5 n and //„ the approximations of the previous section. Fix n > 1. 
Since (p n G L 2 (D), we have n = Ylj a j'^j f° r some eij. Moreover, ||<^ n ||2 = Xlj fl2 an d 
Aoy5 n = — J2j ^j a j¥j- Let Afc( n ) be the eigenvalue nearest to // n . Then 



, 2 a 2 



> (Ajfc( n ) — /i„) 2 a| > (Afe( n ) — /i„) 2 ||^n||2 • 



By (S3D and (IHTLj) . it follows that 



(9.1) |A fc ( n ) - H n \ < 



The right-hand side is a decreasing function of n, so that |Afc( n ) — //J < 0.0987T < ^ whenever 
n > 4. Hence we have the following result. 

Lemma 2. Each interval — j, 1 ^), n > 4, contains an eigenvalue Xk( n )- 
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In particular Xun) are distinct for n > 4. We will now prove that there are only three 
eigenvalues not included in the above lemma. For t > 0, we have (see e.g. [31 [30] ) 

(x, x)dx < I p t (0)dx 



V e -^'= / J2 e ~ X3t (^)) 2 dx= / P?(x,x)dx< / ; 

j=i Jd j=1 Jd Jd 



Tit 



On the other hand, 



N 00 -27rt o 7 

E\ + X "V 717T 4- ts £ I 

e -A fe( „)t > \ e -— * = > 

~ ^ 1 - e~ 2* ~ Tit 2 

n=4 n=4 



for small t > 0. It follows that there are at most 3 eigenvalues of (P/ 3 ) other than \k( n ) 
(n > 4). Furthermore, we have 1 < A x < ^, 2 < A 2 < vr and 3.4 < A 3 < ^ by pQ. 
Therefore, /c(tt,) = n for n > 4, and also by (19. ip . A 3 > 3.83. We have thus proved the 
following theorem. 

Theorem 6. We have 

37T 37T 

1 < Ai < — , 2 < A 2 < vr, 3.83 < A 3 < — , 



2 



and 



n7T 7T 7T ri7T 7T 7T . 

< X n < 1 (n > 4). 

2 8 10 2 8 10 v - ; 



In particular, all eigenvalues of {PP) are simple, \\ n — X m \ > 0.69 when n ^ m, and 

10 



I A n — A TO | > jf: if moreover n > 4. Furthermore, as n — ► oo ; 



7T ~ / 1 

(9.2) A „ = ___ + - 



/717T 7T 



i 

< - n > 1, 



More precisely, 
(9.3) 

v 2 8/ n 

i.e. the constant in O(-) notation in (19.21) is not greater than 1. Indeed, by (19. ip . for- 
mula (19. 3p holds for n > 7, and for n < 6 one can use the estimates (jll.ip . Without 
referring to numerical calculation of upper and lower bounds, one can use (19. ip for n > 4 
and estimates of Ai, A 2 and A 4 of Theorem [6] to obtain (19. 3p with ^ replaced by j-. 
Better numerical bounds for first few eigenvalues are obtained in Section [TTJ 



10. Estimates of eigenfunctions for the interval 

In the preceding two sections the approximations <p n to the eigenfunctions (p n were con- 
structed and it was proved that \i n = ™ — | is close to X n . Now we show that (p n is close to 
ip n in L 2 (D). 

Let n > 4 be fixed. Recall that (p n = Y2j a jfj] with no loss of generality we may assume 
that a n > 0. For j ^ n we have |/i n — Xj\ > j^. Therefore, 

oo 97f^ 
\\A D >p n +Hn<Pr& = ^{^n ~ Xjfa) > (fl n - A n ) 2 a„ + J^^tf- 
j=l j^n 
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Denote the left-hand side by M%; the upper bound for M n is given in (18.31) . We have 

100M„ 2 



Therefore, 

1 1 , _ ... _ _ 20M„ 

V^n - 2 ^ 9 - l^n _ a nVn\\ + { \\<fn 2 _ a n) < 2 ¥>n - <Wn < — 

ii 11/ 

This, together with (18.111) . yields the following result. 



Lemma 3. Let n > 4. With the notation of the previous two sections, we have 1 — < 



fin 



\<p n \\ <l + ±£-, and 



20 / _ 8.00 13.66 1 



\?n-\W»\W\ 2 <—Jl.21 + -—+ . 
In particular, for n > 4, by the above result and (18. lip . 



Wn\ 



20 M n 20 7r _ 2 

3tt ' J 1 _ a^2 < 3tt ' 10 ~ 3' 



Since (p n is symmetric or antisymmetric when n is odd or even respectively, we have the 
alternating type of symmetry of <p n . 

Corollary 3. The function (p n is symmetric when n is odd, and antisymmetric when n is 
even. 

Proof. For n < 3 this is a result of pp. When n > 4, (p n is either symmetric or antisymmetric, 
and the distance between ip n and normed (p n does not exceed |. Therefore (p n has the same 
type of symmetry as (p n . □ 

Corollary 4. As n — > oo ; 

||^-sin((f -|)(l + x) + f)|| 2 = 0^. 



By a rather standard argument, H^nlloc < see e -§- [31] • A slight modification gives 

the following result. 

Proposition 2. Let c = \\<f>n\\o- Then 



C \ V 7T 
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Proof. Let t = Using Cauchy-Schwarz inequality, Plancherel theorem and the inequality 

pf (x, y) <pt(x — y), we obtain 

c\<p n {x)\ <e Xnt \P t D (ap n -<p n )(x) \ +e Xnt \P t D (p n {x)\ 




(p t (x - y)) 2 dy ■ \\c(p n - (p n \\ 2 + y/e \\(p n \ 



e 

2^ 



oo 

e-*WdHI^»-£J + Ve||VwL 

— oo 



e 
2vrt 



— • \\ctp n -ip n \\ + Ve||V/*nlloo' 



and the proposition follows. □ 
Corollary 5. The functions (p n ( x ) are uniformly bounded in n > 1 and x G D. 

More precisely, for n > 1 we have 
(10-2) WL<3. 

Indeed, for n > 7 this follows from (110.11) when the right-hand side is estimated using 
Theorem [6], Lemma [3] and (I5.10p . For n < 6 it is a consequence of Hy^nH^ < y ^ and 

A n S — • 

11. Numerical estimates 

In this section we give numerical estimates for the eigenvalues A n of the semigroup (Pf) when 
D = (—1,1). The following estimates hold true; the upper bounds are given in superscript 
and the lower bounds in subscript: 

A x = 1.157773883697| A 6 = 9.032852690|°|8 

A 2 = 2.75475474221^5 A 7 = 10.6022930996??^ 

(11.1) A 3 = 4.31680106659^| A 8 = 12.1741182627^ 

A 4 = 5.8921474709^ A 9 = 13.744109059^ 

vio = 15.3155549960^ 



A s = 7.460175739^ Am = 15.3155549960P 2 



This is the result of numerical computation of the eigenvalues of 900 x 900 matrices using 
Mathematica 6. 01. Different methods are used for the upper and lower bounds, as is described 
below. For the introduction to the notions of the Green operator and the Green function, 
the reader is referred to e.g. [I]. The explicit formula for the Green function of the interval 
was first obtained by M. Riesz 



11.1. Upper bounds. For the upper bounds, we use the Rayleigh-Ritz method, see e.g. 
Let Gjj be the Green operator for P/ 3 . Then Gr,(p n = j-<p n - The following min-max 
variational characterization of eigenvalues of Gd is well known, see e.g. 



[11.2) — = maxj mini?( / ) : E is //-dnnonsional subspaecof L 2 (D) 
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where R(f) is the Rayleigh quotient for Go, 

j^ 1 f(x)G D f(x)dx 
K {J) = 



Let /„, n — 1,2, be a complete orthonormal system in L 2 (D), and let be the subspace 
spanned by /„, n — 1,2, ...,N. By replacing L 2 (D) by E N in flll.2p . we clearly obtain the 
upper bound X^ N for A„, n = 1, 2, ...,N. On the other hand, (A+jy) -1 is the n-th largest 
eigenvalue of the N x N matrix of coefficients a m n of the operator Go in the basis 
(fi, f2, • /jv) (note that a mjn do not depend on iV). 

The main difficulty is to find a convenient basis /„ for which the approximations converge 
sufficiently fast, while the entries of A N can be computed explicitly. 

For the sake of comparison, recall that analytical computation in [T] gives the upper bound 
Y- ~ 1.178. Our first attempt to use the Rayleigh- Ritz method for Ad instead of G^, with 
fn( x ) — sia( 2 Y(x + 1)), resulted in relatively poor estimates. For example, for iV = 1000 the 
upper bound for the first eigenvalue is A^iooo ~ 1-1579, accurate up to third decimal place. 
A more efficient approach, described below, uses Legendre polynomials. 

We begin with computation the values of the Green operator of the interval (—1, 1) on the 
polynomials g n (x) = x n . Recall that the Green function of the interval D = (—1, 1) for the 
Cauchy process is given by 



, , 1 f i^y? du 1 1 — xy + a/1 — x 2 -v/l — y 2 

GD{x,y) = — — — — = - log ■ 



27T J ^/uy/u+l 71 \X -y 

where x,y G D. Integrating by parts gives, after some simplification, 

G D 9n(y) = J G D (x,y)g n (x)dx 

1 Jl - y 2 f 1 x n+1 dx 
pv ' 



vr n + 1 J_ x s/l - x 2 (x - y) 
1 v/l - y 2 f 1 (x n+1 - y n+1 ) dx 1 y/l - 



7r n + 1 J_i s/l - x 2 (x -y) 7T n + 1 



i(y), 



where 



i{y) = p v j 



dx 



-i Vl - x 2 {x - y) 



The indefinite integral is given by 



log 



\x - 2/1 v 71 - y 2 



xjl^y 2 2(1 - xy + VT^^/T^f 
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and therefore I(y) = 0. Consequently, we have 

, , 1 a/1 -y 2 f 1 (x n+l - y n+l ) dx 
G D 9n{y) ' 



7T 71 + 1 /_! x/i - X 2 (X - y) 

1 v/l - y 2 ^ , afcfe 



2 n /"l 



7T 71 + 1 -f"^ 7-a a/1 - X 2 

n + i r( J + i)- 

Finally, for m, n = 0, 1, 2, ... such that m + n is even we get 

{y)G D g n (y)dy 



— Y , /* a/t^v^- 2 ^ 

" + i^v^r(j + i)J_ 1 v 

i Ul r(j + j) r(^f±i-j) 
2(n + i)f-r(j + i)r(^ + 2-j)- 



By simple induction, one can prove that in this case 

( i r(2^i) r (2±i; 



(11.3) G 



m + n + 2 r(f + i) r(| + i) 



for to, n even, 



i r(f + i) r(f + i) 

for m, n odd. 



[m + n + 2 r(=f2) r(s±2 



If to + 7i is odd, we obviously have G m ^ n = 0. 
The Legendre polynomials are defined by 



LfJ 



n-2i 



fn(x) ^ ^ Cn, 
i=0 

where 

(-l)*(2n-2i)! (-l)T(2n - 2i + 1) 



fll A) c 4 = - 

y n,i 2 n i\(n - i)\(n - 2i)\ 2 n iK(n - i + l)T{n - 2i + 1) ' 

form the orthogonal basis in L 2 (D). Therefore, we have 

-i LfJ Lf J 

(11-5) a m ,n= I fm(y)G D f n (y)dy = Y, Y c n,i c m,jG itj , 

with c nj j and Gjj given by All .3j) and (111.41) . The upper bound for A n is A+^y, where (A^) -1 
is the n-th greatest eigenvalue of the N x N matrix An — (a n>m ). 
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11.2. Lower bounds. To find the lower bounds to the eigenvalues of the problem (11.11) - 
(11. 3p for an interval D = (—1, 1) we apply the Weinstein-Aronszajn method of intermediate 
problems. More precisely, we use (with small changes in the notation) the method described 
in [23] (section The method), where the sloshing problem is considered. For more details, 
see [23] and the references therein. 

The analytic function sin(z) = (sin £ cosh 77, sinh^cos?]), where z = £ + irj, transforms 
the semi-infinite strip R = {(£,77) G R 2 : — f <^<^,?7>0} onto the upper half-space 
H = {(x,y) G R 2 : y > 0}. Let u be a solution to the eigenproblem (ll.ip - (ll.3l) with 
D = (—1,1). Then the image v(z) = u(r)(z)) of the function u under r\ is a solution to the 
following equivalent problem 

(11.6) Av^, V ) = _!<£<! )r? > , 

(11.7) jfr>(£,0) = -Acos£v(f,0) f < f < f , 77 = 0, 

(11.8) v(-*,ri)=v(z i rj) = r]>0. 

For / G L 2 (— ~, |) we denote by A/ (not to be confused with Af) the normal derivative 
of the harmonic function agreeing with / on (— f , f ) and vanishing on {— f,f} x [0, 00) 
(this is an analogue of the Dirichlet-Neumann operator). Since v (£, 77) = sin(A;(£ + ^))e~ kv 



satisfies ( 111.61) and ( 111.81) . the eigenfunctions of A are simply #&(£) = y | sin(A;(£ + |)), and 
We define the operator of multiplication by the function sign £ \/l — cos£ 



(T/)(0 = signal -cos£/(0, / G L 



2^ 7T 7T ^ 
' 2 ' 2 / 



The problem flll.Sp — flll.8p can be written in the operator form as 
(H-9) (Af)(0 = \(l-T 2 )f(0. 

Let Pn be the orthogonal projection of L 2 (D) onto a linear subspace En of L 2 (D) spanned 
by the first iV of the linearly dense set of functions fx, f 2 , .... Then the eigenvalues A"^ of 
the spectral problem 

(11.10) Af = X(l-TP N T)f 

are lower bounds for the eigenvalues of (111.91) and consequently to the eigenvalues X n of the 
problem (|11.6p - (|ll.8p . Roughly, this is because 

7T 

f{x)TP N Tf{x)dx = \\P N Tf{x)\\l < \\Tf{x)\\l = f(x)T 2 f(x)dx, 

and so the Rayleigh quotient associated with (111.101) is dominated by the Rayleigh quotient 
for ( 111.91) . namely 

JT f f(x)Af(x)dx J\ f(x)Af(x)dx 



J_ 2 f f(x)(l - TP N T)f(x)dx f\ f(x)(l - T 2 )f(x)dx 

The problem ( 111.101) is called the intermediate problem. We will later choose /„ so that each 
Tf n is a linear combination of the eigenfunctions of A, say 

K 

(11.11) Tf n = Y,c n ,i9i, n = l,2,...,N, 



i=l 
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where K > N . Let C be the N x K matrix with entries c n> j, and let B be the N x N Gram 
matrix of the functions fx, /jy, i.e. the matrix with entries 

/E. 
2 
I 

Finally, let .D be the K x K diagonal matrix of the first K eigenvalues 1, 2, K of A. Note 
that for each j > K, the function is the solution of (11 1.1 01) with eigenvalue A = j (this 
is because Tgj = 0). On the other hand, if / is the linear combination of gi, g%, Qk with 
coefficients a = (ai, .. then / satisfies (jll.lOp if and only if a is the solution to the 
K x K relative matrix eigenvalue problem, 



11.12) 



Da = X(I -C T B- l C)a. 



By arranging the eigenvalues of (111.12j) and eigenvalues K + l, K+2, ... in the nondecreasing 
order, we obtain the sequence of eigenvalues \~ N of the intermediate problem (111.101) . As 
it was already noted, these are lower bounds for A n . 
We define 

f n (x) = 2y/l + cosxg n (x). 

It follows that 

Tf n (x) = 2smxg n (x) = (-l^g^x) + (-l) n+1 g n+1 (x), 
using the convention that ^o(^) — 0. Consequently, C is iV x (N + 1) matrix of the form 



c 



/° 


1 





■• 











\ 


1 





-1 


■• 
















-1 





1 • • 



















1 


■• 






















•• 





(-1)" 







\o 








•• 


(-1)" 







7 



The coefficients of the Gram matrix B can be easily computed, and we have 



-l) 1+ ^32mn 



ir((m — n) 2 — l)((m + n) 2 — 1) 



+ A5 ri 



whenever m + n is even, and 6 mj „ = otherwise. Finally, the solutions of the spectral 
problem (111.121) are simply the inverses of the eigenvalues of the matrix D~ l (I — C T B~ l C). 



These numbers turn out to be less than N + 2, therefore they form A 



n,JV> 



n 



1,2,.. 



N+l. 



v 



Appendix A. Estimates of p t - p[ 

Let D = (0,oo). Let p t {x,A) = P x {X t e A) for A C R, and fix x > 0. By the strong 
Markov property, 

2P x (X t < 0) = 2P x (X t e D c ) = E x (2p t ^ TD (X(r D ),D c ) ; r D < t). 
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Since 2p s (y, D c ) > 1 for y < 0, s > 0, the right-hand side is bounded below by P x (r D < t). 
Therefore, for t > and x > 0, 

. . . 2 f° t , 2 x ft 

(A.l) P x {td < t) < - / — -. dy = 1 arctan - < mm 1, - 

J -co t 2 + {y-X) 2 7T t \ x 

For t > 0, x,y E D = (0, oo), we have (see formula (2.9) in []]) 
Pt(y~x) -pf(x,y) 1 



t t 



- E x (pt- TD (y - X{t d )) ; T D <t) 



nt Ex {(t-r D r + (y D -X(r D )Y ;TD - t 
< — - P x (t d < t) < mm 



By symmetry, also 

Pt(y-x) -p?(x,y) . f J i_ 

— mm 1 , 

Since — 2) < -7, we conclude that 

/ a r>\ n ^Pt(y-x)-pf(x,y) 1 ./lilt t \ 

(A.2) 0< <-mm , — - , t,x,y>0. 

t 7r \t A x z y z x z y xy z J 

This estimate is used in Sections [2] and [31 

Appendix B. Properties of 77 and B 

In Section [31 a function 77 being the generalized Hilbert transform of — arctan i_ is sought. 
More precisely, rj is the function satisfying 7/(0) = and 

1 /-° 1 



(B.l) n'(t) = - P v / ds, t e R, 

7T J-00 [t - S)(l + S 2 ) 

the integral being the Cauchy principal value when rj < 0. Observe that 

ds = - 1 - / ( -T-—K + — — I 



(t-s)(l + s 2 ) l+t 2 J \l + s 2 t-s 

1 / 1 



t arctan s H — log(l + s ) — log \t — s 



l + t 2 \ 2 

Hence we have 

" (t) = (TT^(2-^ log|i| 

and so 

In particular, 

(B.3) n(-t) = -rj{t) + log Vl + t 2 , t e R. 
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The integrals of over (0, oo) and over (—00, 0) are zero (this follows by a substitution 
u — -), and the maximum and minimum, equal to the Catalan constant C ~ 0.916 and to 
— C respectively, is attained at —1 and 1. It follows that 

(B.4) - log(l + t 2 ) - - < r)(t) < - log(l + t 2 ) + -, te R, 

4 7T 4 7T 

and in particular, 

(B.5) e v{t) ~ y/\F\ as 1*1 °o. 

On the other hand, by (IB. II) . 



1 d r \og\t-s\ A 



7T (it /_„ 1 + S 2 



and for t = 0, 



1 + s 2 Vio Jl / 1 + s Jo 1 + s Jo 1 + s s 



Therefore, 



1 r° iog^ 

irJ-cc 1 + ^ 

A related holomorphic function B plays a major role in Sections [3HZI It is defined by 

1 f log(z- a) 
/ oo 1 + s 5 



(B.6) r)(t) = -l °\ 2 1 ds, t G R. 



(B.7) B{z) = -j °\ J' ds, zeC. 



Here we agree that log(z) = log|z| + i| for 2; G (— 00, 0], i.e. log (and therefore also B) is 
continuous on (—00, 0] when approached from C+, but not from C_; see also Section |H The 
function Re B(z) is harmonic in C \ (—00, 0], continuous in whole C and KeB(t) = r)(t) for 
£ G R. For z e C, we have 



s 



, 1 y iog_^-j ^ 1 r iog(k, , . h 

Re = - y_ - TT ^- d. < - J ^ i + g2 ds = V (\z\), 
and in a similar manner 

1 r° log 1^ — s| ds > 1 r° logHN - s 

7T /_„ 1 + S 2 ~ 7T L 1 + S 



ReB{z) = -l T „ 2 ' ^>- / 1 ' ' 2 1 ds = r/(-|*|)- 



oc 



By (ED, 

(B.8) ^log(l + |^| 2 )-- <R eJ B(^) < Ilog(l + |z| 2 ) + -, zeC. 

4 7T 4 7T 

In particular, 

(B.9) |e B(2) |~v^ as|z|^oo. 

This estimates are used in Section [3] and in Section [7| in contour integration. We also have 
1m B(t) = arctant_; this can be shown directly, or using the first part of this section as 
follows. The function KeB'(t) = i]'(t) is the Hilbert transform of (— arctant_)', and at the 
same time KeB'(t) is the Hilbert transform of — lmB'(t), hence lmB'(t) = (arctan£_)'. 
Since Im_B(0) = = arctanCL, we conclude that lmB(t) = arctant_. 
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The following auxiliary computations related to the functions 77 and B are used in Sec- 
tions [3] and [5j We have 

§ — arctans n 1 2 
as = — arctans (arctans) , 



1 + s 2 2 2 



so that 



1 f 00 % — arctan s , n 

;B10) Hi ^TT^ ds =8- 



By a substitution s 



1 

tan t ' 







- bgd +p ds = _ ^ 

1 + S 2 



We have 



f 2 f 2 f 2 7T log 2 

log sin tar = / log sin tdt + / log cos tdt = / log sin(2t) dt — 

Jo Jo Jo 2 

1 r vrlog2 f 2 . 7rlog2 
— - I log sm wdit — = / log sin udu — . 

2 In 2 In 2 



Therefore, 



Whenever a > — 1 and 6 > ^rp, we have by a substitution 1 + t 2 = - and a formula for the 
beta integral, 

(B.12) / —dt = - / (1 - s)~ s b_ — ds - 2 ; v 



1 + t 2 ) 6 2i v y 2r(6) 
Also, by integration by parts and r(|) = ^/tt, 

poo — tx roo g—tx 

(B.13) / TF —dt = 2x\ -^dt = 2 v ^x~ } x > 0. 

7o t 3/2 Jo Vi 

Appendix C. Estimates for the generator on a piecewise smooth function 
The following estimates are used in Section [HJ Define an auxiliary piecewise C 2 function: 

for x G (—00, 



(C.l) q{x) 



'■> 3/' 



§(x + ±) 2 fori6(-|,0), 
l-f(x-±) 2 forxG(0,§), 
1 for x G (|, 00). 



Note that g(x) + g(— x) = 1. Let / be a piecewise C 2 function on R and let g(x) = q(x)f(x). 
Suppose that g has a compact support. We estimate Ag(x) for x G (—1,0). 

Choose M , M x and M 2 so that |f(x)| < M , |/'(x)| < M x , |/"(x)| < M 2 for x G (-§, §). 
Let 7 = J °° |/(x)|dx. Then 

|/(x)| < M |g"(x)| + 2Mi |g'(x)| + M 2 |g(x)| < 9M + 6M X + M 2 . 
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If z G (— 1, — |), then g(z) = 0, and so Ag(z) is estimated (up to the factor -) by 
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-dx < Mn 



i (a; — z) 

3 v 7 



-rtfe + - 

z 4 



here we used 



g(g) 



i la; — z 

3 v 



< 3M + 



9/ 



^r^)2 < | for x G (— |, |) in the second inequality. For z G (— 1,0) the 
principal value integral in the definition of A can be estimated by splitting it into two parts. 
By Taylor's expansion of g, we have 



pv 



g(x) - g(z) 
(x — z) 2 



dx 



< |-|sup{|/(x)| : xe{z-\,z + \j) 



<|sup{|/(x)| : XG (-|,|)} <3M + 2M! + 2|a; 
for the second inequality note that g"(x) = for x < — g. Furthermore, 



< \9{z 
We conclude that 
(C.2) |^( 2 )| < 

(C.3) |^(^)| < 




dx + 9 \f(x)\dx<6M + 9I. 



3M + |J 



7T 



9M + 2M l + |M 2 + 9/ 



7T 



ZE (-1,0). 
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